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ABSTRACT 

We investigate the collapse of non-spherical substructures, such as sheets and filaments, which are ubiq- 
uitous in molecular clouds. Such non-spherical substructures collapse homologously in their interiors but 
are influenced by an edge effect that causes their edges to be preferentially accelerated. We analytically 
compute the homologous collapse timescales of the interiors of uniform-density, self-gravitating filaments 
and find that the homologous collapse timescale scales linearly with the aspect ratio. The characteristic 
timescale for an edge-driven collapse mode in a filament, however, is shown to have a square-root depen- 
dence on the aspect ratio. For both filaments and circular sheets, we find that selective edge acceleration 
becomes more important with increasing aspect ratio. In general, we find that lower dimensional objects 
and objects with larger aspect ratios have longer collapse timescales. Wc show that estimates for star 
formation rates, based upon gas densities, can be overestimated by an order of magnitude if the geometry 
of a cloud is not taken into account. 

Subject headings: ISM: clouds - ISM: structure - stars: formation 



1. INTRODUCTION 

Molecular clouds are observed to have complex ge- 
ometries and contain non-spheric al su bstructures, includ- 
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ICurrvll2000l: lMversll2009t ). but the collapse properties of 
finite, non-spherical structures have only been considered 
recently and are much less understood, owing to their in- 
herent global and local instabilities, as w ell as the crit- 
ical importance of initial conditions (e.g ., BastienI 19831: 
Bastion et al. 1991 1 iBurkert fc Hartmannil2004l: iHsu et al l 
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observations, in particular, have revealed a plethora of 
filamentary struct ures within star- forming regions (e.g., 
lAndre et al.l l2010l ). Filamentary structures are seen on 
both molecular cloud scales as well as o n the small scales of 
individual protostellar envelopes (e.g., Andre et al.l l2010t 



iTobin et al.l I2OIOI: iHacar fc TafaIIall201lD . Such filamen- 
tary structures are also commonly predicted by star for- 
mation models and formed in molecular cloud simula- 
tions, although their formation mechanisms vary depend- 
ing upon which model is used. Supersonic hydrodynamic 
and magnetohydrodynamic turbulence, gravo-turbulent 
models, gravitational amplification of anisotropics, fiow 
collisions, and global gravitational accelerations have 
all been shown to be capab le of forming filamen- 



Simulations show that non-spherical structures collapse 
on longe r timescales th an equal-density spherica l objects 
(e.g., iBurkert &: Hartman n 2004: Vazquez-Semadeni et al.l 
|2007() and that non-spherical structures are prone to grav- 
itational focusing, whereby strong density enhancements 
form at the edges of these objects (e.g.. iBastienI 
Burkert fc Hartmann 2004; Hartmann fc Burkerd 
Heitsch et al, .2008.: .Hsu et al, 2010.). The longer collapse 
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timescales of non-spherical objects, as well as the pres- 
ence of selective edge accel eration, ha ve also been analyti- 
cally demonstra ted (Burker t fc Hartm ann 2004; P on et ah! 
[2011l:lToala eTal.ll2012l) . 

W hile uniform-density spher es collapse homologously 
(e.g., iBinnev fc Tremain3ll987D , an edge effect in circu- 
lar sheets and filaments produces a second mode through 
which these lower dimensional clouds can collapse, wherein 
the collapse is controlled by an infalling edge sweeping up 
material. It has not been clear, however, whether the edge 
of a lower dimensional structure will obtain sufficient mo- 
mentum to be able to sweep up the interior, so that the 
collapse timescale is controlled by the selective edge accel- 
eration, or whether the roughly homologously collapsing 
interior will slow down the collapsing edge enough such 
that the collapse timescale will approach the homologous 
collapse timescale of the interior. 
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The collapse modes of spherical and infinite, non- 
spherical structures h ave been thoroughly investigated 
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Recently, iToala et al.l ()2012[ ) showed that the free-fall 
times of circular sheet- like ("2D" collapse) and filamen- 
tary clouds ("ID" collapse) depend strongly on the geom- 
etry of the cloud and, for both cases, are larger than that 
of a uniform sphere with the same volume density by a 
factor proportional to the square root of the initial aspect 
ratio, A. For circular sheets, the aspect ratio is given by 
A — R{0)/H, i?(0) being a sheet's initial radius and H 
its thickness, and for filaments, the aspect ratio is defined 
as A — Z{0)/TZ, where Z{0) is a filament's i nitial half- 
length and TZ is the radius. On the other hand, iPon et all 
(|201lh find that the collapse timescale of an infinitely thin 
filament varies linearly with the aspect ratio, rather than 
w ith the square root of the aspect ratio. 

iToala et al.l ()2012f ) obtain their analytic expressions for 
the free-fall times by integrating, over time, the accelera- 
tions at the edges of a circular sheet and cyl inder under 
two d ifferent assumptions. In the first case, iToala et al.l 
(|2012[ ) assume that the mass of the cloud remains con- 
stant and that the density of the cloud remains spatially 
uniform over the entire collapse. This is equivalent to as- 
sumi ng that the c ollaps e is homologous. In their second 
case, IToala et al.l (|2012D assume that the density of the 
cloud remains spatially and temporally constant. For this 
second case, it is assumed that the collapse is dominated 
by the ins weeping edge. These approximati ons are used 
to ma ke the problem analytically tractable. iToala et al.l 
argue that these two approximations represent the 
extreme cases, where the collapse timescale is determined 
solely from either the homologously collapsing interior or 
from the selective edge acceleration, such that the actual 
collapse timescale should lie somewhere between the two 
derived timescales. 

In Section 12.11 we calculate the collapse timescales of a 
uniform-density cylinder, under the assumption that the 
cylinder collapses homologously along its major axis. We 
use the first-order approximation to the acceleration that 
is valid in the interior of the cylinder , rather than usin g 
the acceleration at the edge, as done bv lToala et al.l (|2012f ). 
In Section 1221 we calculate the collapse timescale for a fil- 
ament, under the approximation that the interior density 
remains spatially and temporally constant, by treating the 
gravitational force per unit mass on the edge as a rate of 
momentum transfer p er unit mass , rather than as an ac- 
celeration, as done bv iToala et al.l (|2012[ ). We thus take 
into account the effect of low velocity mass being accu- 
mulated by the edge. We compare our collapse timescales 
to previously obtained results in Section 12.31 In Section 
|3l we compare the relative importance of the homologous 
collapse mode and edge-driven collapse mode in filaments 
and circular sheets. We also discuss in Section |3| whether 
the collapse timescales of filaments and circular sheets have 
the same depende nce upon the aspect ratio, as found by 
iToala et al.l (j2012t ). as well as discussing the implications 
of our collapse timescales. Finally, we summarize our find- 
ings in Sectional 

2. COLLAPSE TIMESCALES OF UNIFORM DENSITY 
CYLINDERS 



While circular sheets and cylinders are formally the 
same type of object, we differentiate between the two based 
upon which axis is longer and along which axis collapse is 
occurring. We refer to objects with radii larger than their 
heights, and collapsing radially, as finite, circular sheets 
and refer to objects with heights larger than their radii, 
and collapsing along their long, vertical axis, as cylinders. 
Thus, finite circular sheets represent "2D" sheets while 
cylinders represent "ID" filaments. 

2.1. Homologous Collapse 

The free-fall collap se timescale of a sphere i s a well- 
studied problem (e.g.. iBinnev &: Tremaindll987() and the 
homologous collapse timescale of the interior of a finite 
circular sheet is calculated by J. A. Toala et al. (2012, in 
preparation [erratum]). For reference, we present deriva- 
tions of the collapse timescales of these two objects in Ap- 
pendices El and |BJ 

We examine a cylinder with a total mass M, a total 
length along the major axis of 2Z{t), and a volume density 
p(t) = M/[2TTTl^Z{t)], where TZ is the time-independent 
radius. We denote the distance of a mass element along 
the major axis from the center of the cylinder as z, the 
initial length of the cylinder as 2Z{0), and the initial vol- 
ume density as p(0). The aspect ratio of the cylinder is 
defined as A — Z{Q)/Tl. The magnitude of the accelera- 
tion along the major axi s of the cylinder, for points w ithin 
the cylinder, is given bv lBurkert fc HartmannI (|2004( ) as 

a{z,t) ^ 2TrGp(t){2z + y^TZ^ + {Z{t) ~ zf 



-y^^T(z(0T^}. ^ ^ (1) 

In a homologous collapse, the density remains spatially 
uniform, although not temporally constant, and all regions 
have the same collapse timescale. For a uniform-density 
object to collapse homologously, the acceleration across 
the object at any given time must be a linear function of 
the radial distance to the collapse center. Equation ([T|) is 
not a linear function of z and thus, the collap s e of a fila- 
ment is not homologous. iBurkert fc HartmannI (|2004l find 
that the first-order approximation to the acceleratioqj, un- 
der the condition that \Z(t) — : 

irGTZ^pit) 



a{z, t) 



> 7^, is 
2z 



z{ty 



(2) 



While Equation ([2|) is also not a linear function of z, for the 
interior of a filament, where z'^ <C Z{t)'^, the acceleration 
can be approximated by 

r 2z 

a{z,t) w 7^G7^^p(^) 



(3) 



Equation (11) of IBurkert fc HartmannI |200l) gives anot her form 
from the first term in the brackets. In our notation, IBurkert~ 

7zGp{t)R^ [- (Zit) + 2)-' + {Z{t) - . 



Since Equation ([3]) is a linear function of z, the assump- 
tion of homologous collapse is reasonable for the interior 
portion of a cylinder. As shown in Appendix [B] this is 
also the case for the interior of a circular sheet. While 
strictly Equation ([3|) should only be used to find acceler- 
ations of the interior of a cylinder, the collapse timescales 
in a homologous collapse are constant across the entire ob- 
ject and thus, evaluating this acceleration at the edge will 
yield the collapse timescale of the interior. In this paper, 
the ends of the major axis of a cylinder are referred to as 
the edges of the cylinder. 

of this first -or der a pproximation, but is missing a negative sign 
li Hartmannl 's I I2004I) Equation (11) should read as |a(z,t)| 
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Equation ([3]) can be re-written in terms of the total mass 
of a cylinder. Evaluating this new expression at the edge, 
where z = Z(t)^ yields 

dv^(i) ^ GM_ 
dt Z{tY' ^ ' 

where wo(0 is the velocity at the edge at time t. This 
differential equation has the same dependences on mass 
and length as in the spherical and circular sheet cases pre- 
sented in Appendices 1X1 and IB] Equation Q can be solved 
for the cylinder collapse timescale, tid, 



(5) 



where tsd is the classical free-fall timescale of a uniform- 
density sphere with the same volume density as the cylin- 
der, as derived in Appendix \K\ Thus, the homologous col- 
lapse timescale for a cylinder is linearly proportional to 
the aspect ratio. 

2.2. Cylindrical Edge Collapse: A Constant Density 
Approximation 

The exact solution for the accelerations of a cylinder. 
Equation ([1]), shows that nonlinear terms become signifi- 
cant toward edges, such that the edges are preferentially 
given more momentum than would be expected for homol- 
ogous collapse. If edges are given sufficient momentum, 
the collapse of a cylinder may be dominated by the edges 
sweeping up interior material, such that the collapse oc- 
curs on a timescale faster than the homologous collapse 
timescale calculated in Section [^?T] 

We now calculate the collapse timescale of a cylinder 
under the approximation that the material inside of the 
edge stays at a constant density and does not move in- 
ward until contacted by the insweeping edge. We assume 
that the only acceleration in the system is due to the grav- 
itational force on the edge caused by the interior material 
and we assume that the edge sweeps up all material it 
contacts, such that the mass of the edge grows as it falls 
inward. We do not consider the gravitational force on one 
edge of the filament due to the other edge. For a cylin- 
der with an aspect ratio of 10, the gravitational force due 
to the second edge only becomes equal to the force due to 
the uniform-density, interior material once the cylinder has 
shrunk to roughly one tenth of its original size. Thus, the 
collapse timescale is relatively unaffected by the presence 
of a mass concent ration at the other edge. As argued by 
iToala et al.l (j2012l) . the true collapse timescale should lie 
between the homologous collapse approximation and this 
constant densit y approximation 

As given by |Pon et al.l (|20lH ) and iToala et al.l ()2012f ), 
the gravitational force per unit mass on the major axis 
and at the edge of a cylinder of length 2Z, radius TZ, and 
uniform-density p is 

(6) 



= 27rGp 2Z{t)+n- ^/W+^Z{rf 



Under the assumption that Z{t) ^ 7?., the square root 
can be expanded and the force per unit mass, to lowest 
order, becomes 

g^2TTGpn. (7) 

Note that the force per unit mass is independent of the 
length of the cylinder. This approximation to the accel- 
eration at the end of a cylinder is within 5% of the exact 



acceleration given by Equation (jS]) for aspect ratios above 
5 and within 10% for an aspect ratio of 3. 

By equating the gravitational force acting on the edge 
with the total rate of change of momentum of the edge, 
we find, in Appendix [Cl that the length of the cylinder at 
a time t is given by 



Zit) « ZiO) 



9J1 
6 



(8) 



where Z{Q) is the initial length of the cylinder. 

Equation ([5]) shows that the effect of the edge of a cylin- 
der accreting low velocity mass is only to lower the effec- 
tive acceleration by a factor of three below the acceleration 
that would be obtained by directly equating the gravita- 
tional force per unit mass to the acceleration of the edge. 
Substituting in the lowest order approximation to the force 
per unit mass at the edge yields a collapse timescale of 



TlD = 



TlD 



' 6Z(0) 
27g7^' 



/32A 



(9) 



(10) 



where, as before, A is the original aspect ratio and t^j^ is 
the classical free-fall timescale of a uniform-density sphere 
with the same volume density as the cylinder. 

2.3. Comparison to Previous Works 

iPon et all (1201 iD examine the accelerations of an in- 
finitely thin, but finitely long, filament and find that the 
acceleration becomes infinite at the edge of such a fil- 
ament . Because of this infinite acceleration, iPon et all 
(|201lh derive a first-order approximation to the acceler- 
ations of their infinitely thin filament that is identical to 
the first-order approximation to the accelerations of a fi- 
nite radius cylinder derived in Section [2] as Equation (j3|). 

We find that the homologous collapse timescale for a 
cylinder is ~ 0.82^1 r^ n. This differs si gnificantly from 
the timescale found bv iToala et al.l ()2012f ) for the homol- 
ogous collapse of a cylinder, ~ 0.92a/At3d- We find that 
the collapse timescale for a homologously collaps ing cylin- 
der scales linearly with the aspect ratio, while iToala et all 
(|2012l ) find a \/~A relationship. The cause of the differ- 
ence between these two results is that iToala et al.l (|2012| ) 
use the acceleration at the edge, while we use a first-order 
approximation to the acceleration in the interior of a cylin- 
der. Thus, we pro be the collapse time scales of the interior 
of a cylinder while iToala et al.l ()2012t ) are sensitive to the 
collapse dynamics at the edge of a cylinder. 

We find that if the collapse of a cylinder is dominated 
by the preferential edge acceleration, such that the inte- 
rior remains static and at a constant density until swept 
up by the edge, t he collapse timescale is yJi2ApK^ t^Y). 
iToala et al.l (|2012[ ) find that the collapse timescale of a 
cylinder, under this constant density approximation, is 
^"ilAj (377^) T3D. Since we show in Section [2?2] that the 
effect of accounting for the additional mass being accreted 
by the edge of the cylinder is to decrease the effective ac- 
celeration by a fac tor of three, it is unsurprising that the 
iToala et al.l (j2012l) collapse timescale is exactly a factor of 
\f7> smaller than what we derive. 
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3. DISCUSSION 

Figure [T] shows the exact accelerations, from Equation 
([l}, of uniform-density cyhnders with various aspect ra- 
tios. The aspect ratios shown range from 2.5 to 20 by 
factors of two. The units of acceleration in Figure [1] are 
1\G jZ ^ where A = irTZ^p, such that the linear approxima- 
tions to the accelerations of all of the cylinders are given 
by the same solid line. The momentum deposited per unit 
length, per unit time of a cylinder can be found by mul- 
tiplying the acceleration by a constant factor and thus, is 
not shown in Figure [T] 

The exact accelerations of a radially collapsing, in- 
finitely thin, unifor-density, circular sheet with mass M, 
radius R(t), and surface density are given by Equa- 

tion (jBl[) in Appendix [B] while the first-order approxima- 
tions to these accelerations are given by Equation (IB2p . 
Figure [5] shows these exact and first-order accelerations of 
a circular sheet in units of AGT,. Figure [2] also shows the 
total momentum deposited per unit angle, per unit radial 
length, per unit time of a circular sheet, as a function of 
enclosed mass, for these two different acceleration equa- 
tions, in units of AGT,R. The momentum deposited per 
unit angle, per unit radial length, per unit time, is related 
to the acceleration via 

p{r,t) = a{r,t)r. (11) 
Figures [1] and [5] show that, for the interiors of uniform- 
density cylinders and circular sheets, the accelerations can 
be well approximated by linear functions of distance. A 
comparison of Equations ^ and (|Bip to Equations ^ 
and (|B2I) shows that this linear approximation is valid to 
within a factor of two for z < 0.7Z in a uniform-density 
cylinder with an aspect ratio of 10 and for r < 0.94i? in a 
uniform-density circular sheet. Near the edges of cylinders 
and circular sheets, nonlinear terms become dominant and 
the accelerations become significantly larger than the first- 
order approximations. This suggests that the collapse of 
cylinders and circular sheets can be described as the com- 
bination of two separate collapse modes. The interiors 
of these objects collapse roughly homologously while the 
outer regions are dominated by an edge-driven collapse 
mode, wherein momentum is preferentially deposited to- 
ward the edges of the objects and collapse proceeds as the 
edges sweep up slower moving interior gas. This prefer- 
ential edge acceleration has been p reviously identified in 

Burkert &: HartmannI 



simul ations and analytic work (e.g.. iBurkert &: Ha 
l2004HHartmann fc Burkertll2007l:1Pon et al.ll20H'l ) 



We associate the collapse timescales of Sections 12.11 and 
I with the two different collapse modes. We consider 
the homologous collapse timescale derived in Section 12.11 
to be the characteristic collapse timescale of the interior 
of a cylinder and the constant density collapse timescale 
derived in Section 12.21 to be the characteristic collapse 
timescale for the edge-driven collapse mode in a cylinder. 

The different dependences on aspect ratio, for a cylin- 
der's collapse timescale, as derived in Sections 12. II and [2.2[ 
is thus explainable as being due to the different calcula- 
tions probing different collapse modes. Our results suggest 
that the interior of a cylinder will collapse on a timescale 
proportional to A while the edge-driven collapse timescale 
will depend on y/A. Furthermore, this suggests that the 
relative importance of the preferential edge acceleration 
will also depend upon the aspect ratio, with the edge col- 



lapse mode being more important in cylinders with larger 
aspect ratios. Such a trend is clearly seen in Figure [TJ 

3.1. Relative Importance of Preferential Edge 
Acceleration versus Homologous Collapse 

While the actual collapse of a cylinder or a circular sheet 
will be a combination of the homologous collapse mode 
and the edge-driven collapse mode, it is possible that one 
of these modes will be dominant. If the edge-driven col- 
lapse mode is dominant, such that the majority of the 
total momentum is injected at the edge, the collapse will 
proceed primarily by the edge falling in and sweeping up 
material. Alternatively, if the majority of the momentum 
imparted to an object is due to the linear term of the accel- 
eration, the infalling edge will have insufficient momentum 
to significantly accelerate the material it sweeps up and the 
collapse will proceed roughly homologously, albeit with a 
slight density enhancement at the edge. 

The total rate of momentum imparted to a circular sheet 
can be calculated by integrating the product of Equation 
(jBip and the mass element over the entire sheet or by mul- 
tiplying the total area under the dashed curve in Figure [H 
by 27r. The total momentum imparted to a cylinder can 
be found by integrating Equation ([T]) over the entire cylin- 
der and then multiplying by the total mass of the cylinder, 
or by multiplying the area under the appropriate curve in 
Figure [T] by the mass of the cylinder. 

For an infinitely thin circular sheet, the total rate of 
momentum imparted to the sheet is finite, whereas for an 
infinitely thin filament, the total rate of momentum im- 
parted is infinite. Thus, there is a maximum limit to the 
fractional contribution of the nonlinear components of the 
acceleration to the total momentum imparted to a circular 
sheet. For an infinitely thin circular sheet, almost twice 
as much momentum is imparted to the sheet due to the 
linear term in the acceleration than due to the nonlinear 
terms. Since the importance of preferential edge accelera- 
tion shrinks with decreasing aspect ratio, the linear accel- 
eration term will dominate the momentum imparted for 
all circular sheets with finite height. Thus, it is expected 
t hat circular sheets shou l d colla pse roughly homologously. 

IBurkert fc HartmannI ()2004l ) simulate the collapse of 
uniform-density circular sheets and find that the collapse 
timescales of the sheets are only 20% shorter than the ho- 
mologous collapse timescale that we calculate. This short- 
ening of the timescale by 20% is almost exactly what would 
be expected given that nonlinear terms contribute an ad- 
ditional 50% to the momentum of in finitely thin circular 
sheets. IBurkert fc HartmannI ()2004t ) also note that the 
interiors of their simulated circular sheets undergo signif- 
icant collapse before encountering the edge, as would be 
expected only if the homologous collapse mode is signifi- 
cant in these circular sheets. 

For uniform-density cylinders with aspect ratios larger 
than five, the nonlinear components of the acceleration 
impart more momentum to the cylinders than the linear 
component. The nonlinear components contribute rela- 
tively more momentum as the aspect ratio increases and 
by an aspect ratio of 10, the nonlinear components con- 
tribute almost twice as much momentum as the linear com- 
po nent. Since realisti c filaments have aspect ratios up to 
60 (jAndre et al.ll20ldl) , preferential edge acceleration may 
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Fig. 1. — Radial accelerations, in units of 2\G/Z, for uniform-density cylinders with different aspect ratios. The dotted, dashed, dash- 
dotted, and dash-triple-dotted lines show the exact radial accelerations for cylinders with aspect ratios of 2.5, 5, 10, and 20, respectively. The 
solid line shows the first-order approximation to the accelerations for all four cylinders. Note how the assumption of homologous collapse 
becomes worse as the aspect ratio increases. 



Froctlonol Enclosed Moss 




Froctionol Rodius 



Fig. 2. — Exact radial accelerations, in units of 4GS, for an infinitely thin, uniform-density, circular sheet are shown as the solid line. The 
first-order approximation to these accelerations is shown as the dashed line. The dotted and dash-dotted lines show the momentum imparted 
per unit angle, per unit radial length, per unit time, in units of AGT?R, based upon the exact and approximate accelerations, respectively. 
The lower x-axis shows the distance along the sheet as a fraction of the total radius and the upper X-axis gives the fraction of the total mass 
of the sheet within that radius. 



6 



control the evolution of many observed filaments. 

3.2. Interpretation 

Equations (|A1I) . (|B4p . and ^ are the differential equa- 
tions describing the collapse of spheres, circular sheets, 
and cylinders, respectively, under the assumption of ho- 
mologous collapse, and all three equations are of the form 

dvo GM 

dt ^ ' ^ ^ 

where x is the collapsing dimension. Because of the 
similarity in these differential equations, the collapse 
timescales of the three different objects are all of the form 



galactic star formation rate, based upon emission from 
and number counts of young, massive stars, yield star for- 
mation r ates close r to a few solar masses per year ( e.g. , 
ISmith et al. 1978: .Diehl et al...2006.: .Misiriotis et al... 20061 



ir^X^ 



When written in terms of the total mass of the cloud, the 
collapse timescales for all three objects are independent of 
the initial aspect ratio. This lack of dependence upon the 
aspect ratio is not trivial, as dimensional arguments place 
no constraints on the proportionality of the unitless aspect 
ratio. 

The dependence of the homologous collapse timescales, 
of spheres, circular sheets, and cylinders, on the aspect ra- 
tio, when written in terms of the initial density, comes 
solely from the conversion between total mass and ini- 
tial density. Thus, the proportionality of the aspect ra- 
tio naturally changes with the changing dimensionality of 
the initial cloud. Each reduction of dimension produces 
an additional y/A dependence in the collapse timescale. 
Since, by definition, A > 1, the "2D" and "ID" collapse 
timescales are larger than that of a sphere with the same 
volume density. Similarly, objects of the same dimension- 
ality but with larger aspect ratios also take longer to col- 
lapse. The increasing importance of an edge collapse mode 
with increasing aspect ratio, in cylindrical structures, will 
partially reduce the difference in collapse times between 
cylinders with different aspect ratios, as well as the dif- 
ference between circular sheets and cylinders. Such an 
edge-driven collapse, however, will still occur on timescales 
longer than the corresponding spherical collapse timescale 
due to the a/A dependence that the edge-driven collapse 
timescale has. 

3.3. Implications 

As discussed by iToala et al.l (|2012|) . the spherical free- 
fall timescale, T3D, is often used to calculate collapse 
timescales, and thus star formation rates, from observed 
gas densities, regardless of geometry. Since recent obser- 
vations reveal a multitude of non-spherical substructures 
withi n molecular clouds ( e.g., lMver£|[2009t iMolinari et al.l 
l2010t lAndre et all |2010[ ). using the spherical free-fall 
timescale underestimates collapse timescales and overes- 
timates star formation rates. In particular, attempts 
to predict the total galactic star formation rate from 
the observed gas properties of the Milky Way, by di- 
viding the total molecular mass of the galactic interstel- 
lar medium by the spherical free-fall time correspond- 
ing to the mean density and temperature of the molec- 
ular gas, produce va l ues o f at least 2>Q Mp, ^^'^ 
IZuckerman fc Palmeil 11974 IZuckerman fc Evanj 11974 
while more direct, observational determinations of the 



Murrav fc Rahnranll2010l: iRobitaille fc Whitnevll2oTo[ r 

For a filamentary structure with an aspect ratio of 60, 
correspondin g to the upper lim it of observed filamentary 
aspect ratios ()Andre et aLl|20 ltf), the homologous collapse 
timescale is a factor of almost 50 slower than the cor- 
responding spherical free-fall time. The edge-driven col- 
lapse timescale, while faster than the homologous collapse 
timescale, is still almost 15 times slower than the spher- 
ical free-fall time. For a more typical aspect ratio of 10, 
both the homologous and edge-driven collapse timescales 
of such a filament would be roughly seven times slower 
than the spherical free-fall timescale. Thus, geometric 
considerations can account for a considerable portion of 
the discrepancy between observed and predicted galactic 
star formation rates, although they are unlikely to ac- 
count for the entire discrepancy. For a further discussion 
on the implications of "ID" and "2 D" c ollapsing objects 
having long e r time scales, please see iPon et al. (2011.) and 
iToala et al.l (|2012[ ). 

The two collapse modes studied in this paper, the ho- 
mologous and edge-driven collapse modes, are both global 
collapse modes. That is, these modes cause a cloud to col- 
lapse into one central object. Density perturbations within 
molecular clouds will introduce local collapse modes and 
these local collapse modes must operate on timescales less 
than the global collapse modes, as molecular clouds are 
observed to fragment into clusters of stars, rather than 
c ollapsing to form million solar mass stars. 

IPon et all (l20Tll) examine the conditions under which 
local collapse modes are significantly faster than the ho- 
mologous global collapse mode in spheres, circular sheets, 
and cylinders. They find that strong perturbations are re- 
quired for local collapse modes to be significantly faster 
in circular sheets and spheres, but small (~10%) den- 
sity perturbations in a cylinder can collapse significantly 
(three times) faster than the entire cylinder if the total 
length of the cylinder is greater than 10 times the length 
of the perturbation. Since thermal motions support per- 
turbations smaller than the Jeans length and the radial 
length scale of a radially supported cylinde r is also approx- 
imate l y the Jean s lengt h (Stod61kiewic'3 Il963t lOstrikeij 
I1964D . IPon et al.l (|201lD find that local collapse modes 
are most effective in high aspect ratio cylinders. Unfor- 
tunately, it is in these large aspect ratio cylinders that 
the preferential edge acceleration produces significantly 
faster collapse timescales than the homologous collapse 
timescale. Thus, cylinders, as well as circular sheets and 
spheres, may require strong density perturbations or large- 
scale support mechanisms for local collapse modes to be 
significantly faster than global collapse modes. Note, how- 
ever, that preferential edge acceleration naturally produces 
strong density perturbations along the edges of circular 
sheets and fila ments, as seen in simulations of both circu- 
lar s heets (e.g.. iBurkert fc Hartmannii2004i ) and filaments 
(e.g., lBastienlil983ir ^ 



3.4. Caveats 



Aspect Ratio Dependence of Collapse Times 



7 



Realistic cloud structures are not perfect uniform- 
density spheres, circular sheets, nor filaments, and devi- 
ations from such perfect, symmetrical shapes will influ- 
ence the collapse pro p erties of the clouds. For instance, 
iBurkert fc HartmannI (|2004D show, via simulations, that 
deviations from axisymmetry produce gravitational fo- 
cusing, whereby local density enhancements are readily 
formed. 

For the timescales calculated in this paper, it is assumed 
that clouds have sharp density boundaries. More real- 
istic clouds are likely to taper off slowly at the edges. 
The introduction of such density tapers at the edges of 
clouds is known to reduce the significance of the prefer- 
ential edge acceleration, although such a taper has to be 
quite large in comparison to the size of the constant den- 
sity interior before the preferential edge acceleration is sig- 
nificantly weake ned (Nelson & Papaloizou 1993; Li 2001; 
iPon et aLll201l[ ). As such, preferential edge acceleration 
may be slightly weaker in realistic clouds than assumed 
here, but edge-driven collapse modes should still be im- 
portant for reasonably elongated filaments. 

Star-forming regions are observed to be turbulent 
and this effect can pr ovide support against grav- 
itational collapse ( e.g., iHennebelle fc Chabrieil l201lt 
iPadoan fc Nordliindl |2011|) . Thermal pressure, rotation, 
and magnetic fields can also significantly alter the col- 
lapse of a cloud. Thermal pressure, however, is only ef- 
fective at supporting objects on scales smaller than the 
Jeans length, whereas observed filam ents are often muc h 
longer than the Jeans leng th (e.g., lAndre et al.l |2010|) . 
iBurkert &: HartmannI (|2004l ) also point out that, if solid 
body rotation were supporting the interior of non-spherical 
clouds, the edge would still collapse due to the preferen- 
tial edge acceleration, and if the exterior were rotation- 
ally supported, the central regions would be moving too 
rapidly and would expand. Finally, magnetic fields are 
only capable of providing support perpendicular to the 
field lines and thus, may not be capable of supporting fila- 
ments and circular sheets depending upon the orientation 
of the magnetic field. Models of filament formation predict 
that magnetic fields can be either parallel or perpen dicu- 
lar to the long axis of a filament (jNagai et al.lll998l ) and 
observational studies have found magnetic fields that are 
both perpendicular a nd parallel to the long axes of fil- 
amen t ary structures (iGoodman et al.l 19901 Houde et al 



2004 IVallee fc Fieed 120061 l2007t iSchneider et al 
Chapman et al.ll201ll: ISugitani et al.ll201l[ ) 



et al.l 

2oToF 



The influence of collapse along the short axis of any 
cloud has not been considered because thermal support 
will generally be more effective at preventing collapse along 
shorter axes. In deriving the collapse timescales, the as- 
sumption that A 3> 1 has been utilized. While this as- 
sumption clearly breaks down at late times in the collapse, 
the collapse timescale should be primarily dependent upon 
the early stages of the collapse when infall velocities are 
still relatively small. Local collapse modes are also ex- 
pected at later times in the collapse when, as a conse- 
quence of the density increase caused by the global col- 
lapse, the short dimension becomes larger than the Jeans 
length. 

The cases considered in this paper, pure homologous col- 
lapse and pure edge-driven collapse, are idealized collapse 



modes and are likely to bracket the true mode of collapse. 

4. SUMMARY AND CONCLUSIONS 

We have calculated homologous collapse timescales for 
the interiors of uniform-density cylinders, based upon a 
first-order approximation to the accelerations along the 
major axis of the cylinders. We have also calculated the 
collapse timescale for a uniform-density cylinder, under 
the approximation that the central density remains con- 
stant, by associating the gravitational force per unit mass 
on the edge with the rate of change of the momentum per 
unit mass of the edge. With these results, in conjunc- 
tion with the homologous collapse timescales of uniform- 
density circular sheets calculated by J. A. Toala et al. 
(2012, in preparation [erratum]) and reproduced in Ap- 
pendix jB] we find that 

• Two separate collapse modes are present within 
circular sheets and filaments. The interiors of 
these clouds collapse roughly homologously while 
the edges are preferentially given more momentum, 
such that the edges sweep up material and form 
density enhancements. The effect of preferential 
edge acceleration has been previousl y noted i n sim- 
ulations and analytic studi e s (e.g.. iBastie n 1983; 
Burkert fc HartmannI 2004t Hartmann &: Burkert 



2007 



2008; 



2012) 



Vazquez-Semadeni et al.l 2007t Heitsch et al. 



Hsu et al.ll2010t iPon et al.l 120111: iToala et al 



• The homologous collapse mode is dominant in circu- 
lar sheets while the edge-driven collapse mode dom- 
inates the momentum imparted to filaments with 
aspect ratios larger than 5. 

• The homologous collapse timescales for the interi- 
ors of filamentary clouds scale linearly with A = 
Z{0)/TZ, where 2Z{0) is the total initial length and 
TZ is the radius of the filamentary cloud. The edge- 
driven collapse mode (constant density collapse) 
of filamentary clouds produces collapse timescales 
that are proportional to ^/A. Thus, preferential 
edge acceleration is most important for clouds with 
large aspect ratios. 

• Regardless of dimensionality, the acceleration, un- 
der the assumption of homologous collapse, can be 
expressed as 

dvo GM 

with X ^ the collapsing dimension (e.g., radius R 
for spherical and sheet-like clouds, and the semima- 
jor axis Z for a cylinder). Thus, each reduction of 
dimension produces an additional ^/A dependence 
in the homologous collapse time. 

• In general, lower dimensional objects ("ID" and 
"2D") and objects with higher aspect ratios have 
larger collapse timescales than for a sphere with the 
same volume density. 

• Estimates of star formation rates from gas densities 
can be overestimated by an order of magnitude, for 
realistic filamentary aspect ratios, if the geometry 
of a cloud is not taken into account. 
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APPENDIX 

APPENDIX A: HOMOLOGOUS COLLAPSE TIMESCALE OF A SPHERE 

The free-fall collapse timescale of a uniform-density sphere is a well-studied problem fe.g.. iBinnev fc Tremainelll987^ 
and it is known that the collapse proceeds homologously. For a sphere with a volume density of the acceleration at a 
distance r from the center of the sphere, a{r) = AnGpr/S, is linearly dependent upon the radial distance, as required for 
homologous collapse. 

The governing differential equation for the collapse of a sphere is 

dvo{t) _ CM 

where VQ{t) is the velocity of the edge at time t, R{t) is the total radius of the sphere at time i, and M = 47rp(0)i?(0)"^/3 
is the total mass of the sphere, with p(0) being the volume density at i — 0. This differential equation is well known and 
can be solved for the classical free-fall timescale of a sphere, tsd, 



APPENDIX B: HOMOLOGOUS COLLAPSE TIMESCALE OF A CIRCULAR SHEET 

The homologous collapse timescale of a circular sheet is derived in J. A. Toala et al. (2012, in preparation [erratum]) 
and the general derivation is reproduced here. We consider a circular sheet with mass M, radius R.{t), a nd surface density 
S(t), s uch that the initial radius and surface density are R{0) and S(0) — Af/[7ri?(0)^], respectively. Burkert fc HartmannI 
(|2004[ ) show that for an infinitesimally thin circular sheet, the radial acceleration at a distance r from the center of the 
sheet is 

„,,,, = 4GEWS2[,,(_L_)_E , (Bl) 

where K is the first complete elliptic integral and E is the second complete elliptic integral. They note that the acceleration 
of such an infinitely thin sheet is infinite at the edge, where r — but that in a sheet with finite thickness, the 

acceleration at the edge is finite. 

Du e to this infinite acceleration at the edge, it is common (e.g.. lBurkert fc Hartmannll2Q04l : [Pon et aLll201ll : lToala et all 
I2OI2I) to use a first-order approximation to the acceleration, 

a(r,i)«vrGS(t)-^. (B2) 

The terms excluded from the above equation are higher order terms of r / R{t), such that the above equation is not valid 
at, or near, the edge of a circular sheet, but is reasonably accurate for most points within the interior of a circular sheet. 

While we use Equation (|B2I) in the following derivation, we assume that the sheet has a constant, finite height of H 
that is much smaller than the initial radius i?(0), such that we can define the aspect ratio as A = R{0)/H. We thus 
rewrite Equation (|B2[) as 

air,t)^7rGp{t)Hj^^, (B3) 

where p{t) = Yj(t)/H is the volume density of the sheet at time t. 

It is critical to note that the acceleration in Equation (jB3p is a linear function of the radius and, thus, this equation 
exactly describes a homologous collapse. That is, to first order, the accelerations of an infinitely thin, circular sheet will 
cause the sheet to collapse homologously. It is only the higher order terms oir/R that cause a deviation from homologous 
collapse. 

Denoting the velocity at the edge as VQ{t) and using the relation between the mass of a circular sheet and its volume 
density, M — irR^p/A, Equation (jB3|) can be re-written as 

'^^o ^ GM 

1^~rW ^^^^ 

Equation (|B4p can be solved to show that the collapse timescale of a circular sheet, T2d, is 



T2DWY-^T3D. (B5) 

Thus, under the assumption of homologous collapse, the timescale for collapse of a circular sheet scales with \/~A. 

APPENDIX C: UNIFORM-DENSITY COLLAPSING FILAMENT 

Assuming that the only force on the edge of a cylinder is from the uniform-density interior material, this gravitational 
force from the interior must be equal to the rate of change of the momentum of the edge. From Newton's second law, 

m{t)g ^v{t)—^ +m{t)—^, (CI) 
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where v{t) is the velocity of the edge and m(t) is the mass of the edge. We define the sign of the velocity such that v{t) 
is positive for inward motions. See also Section [2^ for other variable definitions. Since the mass at the edge increases as 
material is swept up 

^-v^7^V.W. (C2) 
Thus, the differential equation governing the motion of the edge is 

m(^)5 = 7r7^2p^;(^)2+m(^)^^. (C3) 

at 

For the remainder of this derivation, we drop the functional dependence on t from our notation. Rewriting the derivative 
of the velocity yields 

dv dv dm 
dt dm dt ' 

— ^ —nn^pv, (C5) 
dt dm 

mg — ttTZ^ pv^ + m— — ttTZ^ pv, (C6) 
dm 

dm m TT/c^ p 



(C4) 



w 

9 ' 



This is an Abel differential equation of the second kind. It can be solved with the aid of the substitution 

(C8) 
(C9) 
(CIO) 

— • (Cll) 

At the beginning of the collapse, m = and v = 0, such that w{t = 0) = 0. We have also used the fact that the 
gravitational force per unit mass, to lowest order, is constant over the collapse. 
The edge velocity, as a function of the edge mass, is thus 







w 




V = 








m ' 




dv 


dw 1 




dm 


dm m 


dw 


w 
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dm 


m^ 


nW^p' 


w{t = 


0)2 _ 


(nr' — mi 



37r7c"^ p 

Since the mass at the edge is equal to the mass swept up. 



m = (Z(0) - Z)7r7^V, (CIS) 
dZ 



dt 



'2(Z(0) - 


z)g 


3 


k{Z{{))- 


z) 



(C14) 



i=y^^^ (C15) 

Z = Z(0)-^. (C16) 



